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r* Abstract. The standard theoretical description of coherent backscattering, accord- 
Ch ing to which maximally crossed diagrams accounting for interference between counter- 
ed propagating path amplitudes are added on top of the incoherent background, violates 

iw the fundamental condition of flux conservation. In contrast to predictions of previous 

, t theories, we show that including maximally crossed diagrams with one additional scat- 
tering event does not restore flux conservation. Instead, we propose that the latter is 

C recovered when treating the effects of coherent backscattering and weak localisation in 

i-£j a unified framework. On the basis of this framework, we demonstrate explicitly flux 

£JJ conservation in leading order of the weak disorder parameter l/(k£). 

o 
o 

PACS numbers: 42.25.Dd, 73.20.Fz, 42.25.Fx 
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Q\ 1. Introduction 

00 

The effect of coherent backscattering leads to an enhancement of intensity in exact 
backscattering direction when measuring the average wave intensity scattered from a 
random medium. Known since the years 1984/85 where it has been both described 
theoretically [U |2] and observed in experiment [3j Hj, the acknowledged theoretical 
explanation commonly presented in today's textbooks (see e.g. [3]) holds the following 
mechanism responsible for this phenomenon: The waves propagating within the medium 
collect random phases due to the random realization of the scattering potential. Thus, 
upon averaging over different realizations of the disorder, only those constellations 
of the intensity propagator survive for which the phase difference between the wave 
and its complex conjugate counterpart vanishes. Two possible kinds of propagation 
processes are known to lead to exact phase cancellation: The case where the wave and its 
complex conjugate visit the same scatterers in the same order (ladder propagation); they 
propagate the same distance at each step, thereby collecting equal phase. With regard 
to backscattering, this process leads to a background intensity which is distributed 
smoothly over all backscattering angles $, since phase cancellation occurs for every 
backscattering angle (incoherent background). However, for systems which exhibit 
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reciprocity symmetry [6], exact phase cancellation is also achieved in the case of counter 
propagating waves (maximally crossed propagation), where again the same scatterers 
are visited by the wave and the complex conjugate, but in reversed order. Resulting 
in equal phases in exact backscattering direction, i.e. backscattering angle i? = 0, but 
leaving a nonzero random phase shift for every angle other than zero, this gives rise 
to the sharp backscattering peak centered around $ = 0, known as the characteristic 
coherent backscattering peak. 

This explanation, however, implies a serious deficiency: it is possible to show [7] 
that for a system where only ladder-like propagation processes are considered, the 
incoming intensity flux (f in of a wave entering a medium is equal to the outgoing 
intensity flux (p ou t, and hence flux conservation is fulfilled. Adding the coherent 
backscattering peak - which yields a non-negative contribution to the outgoing flux for 
every backscattering angle - on top of the incoherent background therefore violates the 
fundamental law of flux conservation. Therefore, we may ask which scattering processes 
counterbalance the maximally crossed contribution and thereby restore flux conservation 
in the process of coherent backscattering. This question has been addressed before |8] , 
where experimental studies are described together with a brief theory to explain the 
results obtained. 

In the work presented here, after briefly outlining the theoretical frame (see 
chapter [2]), we will show that the ansatz proposed in [8j, i.e. including maximally 
crossed scattering processes which contain one additional scattering event, does not 
restore flux conservation (see chapter pi). Then, we will present the mechanism which 
we conjecture to meet this goal instead: We consider the full set of possible loop 
propagation processes by accounting explicitly for the presence of a boundary surface in 
the case of a finite or semi-infinite medium (see chapter El). This amounts to a complete 
and consistent treatment of coherent backscattering together with the effect of weak 
localisation [9] for wave propagation in disordered media. Finally we verify that this 
consistent treatment indeed restores flux conservation for coherent backscattering in 
leading order with respect to the disorder strength l/(k£) (see chapter [5]). 

2. Theoretical frame 

We treat the propagation of a monochromatic, scalar wave through a disordered medium 
where it undergoes multiple random scattering. The disorder is described by a potential 
V^(r) denoting a random function of position r. For the disorder average over many 
realizations of the disorder of any quantity X we write X. We work in the frame of the 
Gaussian white noise model being fully characterized by its zero mean value V(r) = 
and non-zero correlation function V(r)V(r/) = B(r, r/) as well as vanishing cumulants of 
all orders higher than two. Assuming both, rotational as well as translational invariance 
of the correlation function, we may write B(r, r/) = B(\r— r/|). Furthermore, introducing 
£ c as the correlation length of the potential, i. e. as the characteristic length of the decay 
of the correlation function, we restrict ourselves in the following to the regime where 



the wavelength A of the scattered wave is much larger than the disorder correlation 
length, i.e. to the regime A 3> £ c . In this case, the Dirac delta distribution is a good 
approximation for the correlation function for which we hence write 

B(\r-r/\) = u 2 5(\r-rf\) , (1) 

with u 2 as a pre-factor taking the role of the scattering strength of the potential. 

To describe the propagation of a wave with wave number k in the presence of the 
potential V(r) we are in general looking for wave amplitudes ip(r) satisfying the scalar 
Helmholtz equation 

(A + k 2 (l + V(v)))4>(r)=p(v), (2) 

for given source distribution p(r) where A = V 2 denotes the Laplace operator. 

For the vacuum case (i.e. V = 0) the Green's function for the Helmholtz operator 
solving pi) for a delta-like source p(r) = 5(r — r/) is given by the well-known (retarded) 
vacuum Green's function [3] 

ik\r—rt\ 

Gf = G„(r-«) = - sjF - sf (3) 

or by the advanced vacuum Green's function Gq = (Gq)* as the complex conjugate of 
pi). Go(r — r/) describes the spherical wave at point r which emanates from a point-like 
source at r/ having propagated the distance |r — r/| through empty space. 

Likewise it is possible to determine the disorder averaged Green's function G(r — rf) 
to (J2j) in the presence of a disorder potential V(r) ^ 0, where taking the disorder 
average again implies that instead of specifying the solution for one specific realization 
of V(r) we describe propagation between the points r/ and r averaged over many different 
configurations of the disorder. G in spatial representation reads [5] 

p ifc|r-r/| 

and exhibits the same functional form as the vacuum Green's function Go, eq. (J3J), 
except for the wave number k of fe) which is replaced by k in Ml). This effective wave 
number k turns out to be a complex quantity the non-zero imaginary part of which 
leads to an exponential damping of the Green's function as a consequence of scattering 
within the disordered medium. The corresponding decay constant l s is determined by 
the imaginary part of k in the following way: 

h - k+ k (5) 

l s is referred to as the scattering mean free path describing the average distance between 
two successive scattering events. For simplicity, we assume that the real part of k is 
unchanged. (If necessary, this can be achieved by adding a constant term to the potential 
V) 

The scattering mean free path certainly depends on the properties of the disordered 
medium and the scattering processes. In the case of sufficiently dilute media, i.e. if the 
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wavelength of the scattered wave A = 27r/fc is much smaller than the scattering mean free 
path l s , it is justified to assume only the above described ladder propagation processes to 
survive disorder average. This approximation is called ladder approximation and leads 
to the following expression for the scattering mean free path [5]: 

Arc 

ls,Ladder = £ = ~ T- (6) 

Hence, within the frame of ladder approximation, the scattering mean free path depends 
only on the scattering strength u 2 of the disorder potential as it was introduced in Q. 
Corrections of l s beyond the ladder approximation will be considered in chapter [5] 
Furthermore, the expression hi turns out to be a system parameter classifying the 
strength of the disorder by comparing the wavelength of the scattered wave to the 
scattering mean free path. Hence, a large value of the disorder parameter (k£ ^> 1 •<=>- 
£ ^> X/2tt) signifies weak disorder. 

3. Analysis of different scattering scenarios 

Following the idea of |Sj, according to which the contribution of the maximally crossed 
scattering sequences (H ^-contribution in figure [I] below) could be cancelled by the 
contributions of processes which contain one additional scattering event, we analyze in 
this chapter the contributions of all such scattering sequences one may possibly construct 
(see figure [I]) . 

For the analysis we use the bistatic coefficient [10] 7(f2) = <fout(^) / >( fin as the ratio 
of incoming and outgoing intensity flux which is defined as follows: We assume detection 
of the scattered intensity at a distance R from the scattering medium sufficiently large 
to allow far field description and for the emergent wave to be of approximately spherical 
shape. Therefore, we write 

yw(n) = lim /(i?0)47ri? 2 , (7) 

(fin = AI , (8) 



where I(Rft) = \^fj(RQ,)\ 2 denotes the average wave intensity at point R = RQ, and 
Io and A refer to the intensity of the incoming wave and the surface of incidence, 
respectively. Combining eqs. (J7|) and (|8j we find the bistatic coefficient to read 

AttR 2 
7 (0) = lim -jrj-HRn), (9) 

Integrating over all angles f2, the condition for flux conservation thus reads: 

^7(n) = i (io) 

Using this definition, we now calculate the bistatic coefficients for the scattering 
processes shown in Fig. [TJ Each diagram contains incoming and outgoing lines associated 
to the points ri, r 2 or r 3 , described by the functions ipi n and G out , see eqs. ( 12|l3 ) 



below, and propagators connecting these points with each other, which will be discussed 
separately below. 




Figure 1. Scattering processes investigated beyond ladder approximation. Solid lines 
represent Green's functions and dashed lines the respective complex conjugates. The 
case of maximally crossed propagation (Ha) is responsible for the famous coherent 
backscattering cone which, considered separately, violates flux conservation. In order 
to unearth possible cancellation of this crossed contribution, we construct further 
scattering scenarios which contain one additional scattering event at point r3. We 
take into account Hg (or (He)* = He, for complex conjugation exchange every solid 
line by a dashed line and vice versa) and Lp, (or (Lp)* = Lp) where the crossed 
propagation process between points ri and T2 is replaced by a ladder process. Since 



ladder and crossed propagation are equivalent (see (17)), these four diagrams yield 
equal contributions. To arrive at process Lp (or (Lp)* = Lq) we interchange dashed 
and solid lines to construct a mixed event where the additional scatterer is visited once 
by the wave and once by its complex conjugate. L single (or (Lsingie)*) corresponds 
to a degenerate Hb or Hq process where the points ri and V2 are merged, i.e. no 
scattering sequence takes place between these points. The backscattering contributions 
originating from these scattering processes, as described by eqs. (fl8j [20l [22] [24] ), are 
analyzed, first as a function of the backscattering angle $, see figure [3] and second 
integrated over all backscattering angles to obtain the total backscattered flux, see 



eqs. (25-28) 



For the calculations we consider a slab geometry exhibiting finite thickness L in 
^-direction but extended to infinity in x and y-direction (see figure [2]). This implies 
A — > 00 in (J9J), which, however, will be counterbalanced in the following by restricting 
integrations over the whole scattering medium to integrations over z: 

dr^ [ dz, (11) 
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The optical thickness of the medium is given by b = L/L We assume the incoming wave 
to be plane with wave vector k in = ke z perpendicular to the x-y-plane as the surface of 
incidence at z = 0: 



JKin-r 



Akz 



Tj, in (r) = e ~^e lK ^ r = e~^e lKZ , (12) 

The outgoing wave is described by the average Green's function as given in Q in 
Fraunhofer approximation [5]: 



G out (r,R) = e 



3 ,5 e -ik\R.-r\ 

47r|R-r| 



2/ COS # .. g 

W s g — tk-out-r 



ikR 



AllR' 



(13) 



with k oui = kH/R and z/ cosd the distance the wave has to travel from the point 
r until leaving the medium. Reflection at the boundary of the scattering medium is 
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Figure 2. Setup and coordinates chosen for the calculations: We consider a slab 
geometry of finite thickness L in z direction but extended to infinity in the x and 
y direction. The incoming wave is assumed to be a plane wave propagating in z 
direction, thus entering the medium perpendicularly to the surface of incidence. The 
backscattering angle i? is defined as the angle between the incoming and outgoing wave 
vector: # = <(k cmt ,kj„). 



neglected for, both, incoming and outgoing waves (which is appropriate in the case of 
weak disorder). Due to rotational symmetry around the z-axis, the bistatic coefficient 
for the slab geometry only depends on the backscattering angle $, see figure |2j such 
that the condition of flux conservation, see eq. (10), turns into J_ 1 ° s 7(1?) = 1. 



Concerning the remaining propagators in figure [TJ the curly lines connecting the 
points i - ! and r 2 represent the crossed (for diagrams Ha and H B ) or ladder (L F 
and Ld) propagator Pc(^ 1 ,r 2 ) or Px,(ri,r 2 ) for the average intensity. The ladder 
propagator fulfills the following self-consistent integral equation, describing a ladder- 
like propagation process between the points ri and r 2 [2]: 



^L(ri,r 2 ) = Po,z(ri,r 2 ) 



dr i P L (r 1 ,r i )P 0:L (r i ,r 2 ) 



(14) 



where ri stands for the positions of all possible intermediate scattering centers we 
integrate over and Po,l denotes the single-step ladder propagator given by 

e -\n-r 2 \/i 



-Po,i(ri,r 2 ) =u G(r 1 ,r 2 )G (ri,r 2 ) 



(15) 



47r^|ri -r 2 | 2 ' 

with l s = £ in eqs. (4|5). Eq. (14) can be interpreted as a random walk of a classical 
particle with average step length £ and isotropic scattering events between the individual 
steps. The corresponding bistatic coefficient results as 



7 ' Ladder (&) = ^R 2 



c^ m (,2)|V|G out (r,R)| 2 



+ / dz 1 dY 2 \i) in (z 1 )\ 2 P L (r 1 ,Y 2 )u 2 \G out {r 2l K) 



(16) 

(with the first term describing single scattering) and fulfills flux conservation, i.e. 
j dc ° s ^ 'y Ladder ^'j = 1 [7]. Hence, all additional processes added on top of the 
incoherent ladder background ^Ladderi.^) have to cancel each other mutually in order 
to maintain flux conservation. 

The crossed propagators in figure [T] describes a process where the two conjugate 
amplitudes propagate in different directions (one from r x to r 2 , the other one from r 2 to 



ri). However, since the Green's function is symmetric under exchange of its arguments, 
the resulting propagator is equivalent to the ladder propagator 

PL(ri,r 2 ) = P c ( ri ,r 2 ), (17) 

We now consider separately the different scattering processes shown in figure [I] 

Ha- Crossed contribution process which is responsible for the effect of coherent 
backscattering. Considering its contribution, only, entails violation of flux 
conservation. We write for the corresponding bistatic coefficient 

1h a (i?) = 47i\R 2 / / dz 1 dr 2 i/j in (zi)ilj* n (z2)PH A (ri,r2)u 2 G ou t(r2,K)G* out (ri,K) 



\ [ /dz 1 rfr 2 e~ ( ^ 1+ ^ )( ^ + ^^ ) P c (r 1 ,r 2 )e i(k " l+kout)(ri " r2) , 



where we used the explicit forms of G out and ip in given in eqs. (12) and (13), 
respectively (with l s = £), as well as the fact that the propagator for a crossed 
propagation sequence Ph a m &y be identified as 

P^(r,r/)=P c (r,r/). (19) 

Hb- Describing a scattering sequence of crossed propagation of wave and complex 
conjugate now supplemented by an additional scattering event with incoming and 
outgoing wave encounter at position r 3 : 

7h s (^)=47tP 2 / / / dzidr 2 dr 3 ilj in (z 3 )ilj* n (z 1 )PH B (r 1} r2 } r 3 )u 2 G ou t(r 3} K)G* out (r2,K) 



1 / / ^^ 1 rfr 2 rfr 3 e- ifl ^e-^^P^ s (r 1 ,r 2 ,r3)e ikm - (r3 - ri) e ikot ' t{r2 - r3) , 




(20) 

where Ph b denotes the propagator corresponding to this particular scattering 
sequence, 

P^(r l7 r 2 ,r 3 ) = G(r 3 ,r 2 ) M 2 P L (r 2 ,r 1 )G(r 1 ,r 3 ). (21) 

In [8], the authors claim that this process - together with its complex conjugate 
7-HcC^) = ThbW _ restores flux conservation, i.e. that it counterbalances the 
coherent backscattering contribution 7# A ($). 

Lp\ A scattering process containing one single additional scattering event at point r 3 
as in the Hz-case, but now visited once by the wave and once by the complex 
conjugate: 

1l f (&) = 47rP 2 / / dzidr 2 dr 3 tfj in (z 3 )ilj* n (z 1 )P LF (r 1 ,r2,r 3 )u 2 G ou t(r2 J K)G* out (r 3 ,K) 



' (/ / ^^ 1 t/r 2 t/r 3 e" i£ ^e-^^P LF (r 1 ,r 2 ,r 3 )e ik " l - (r3 - ri) e ikout(r3 " r2) , 




(22) 
with P Lf the propagator for this scattering sequence, 

P L A r u r 2,r 3 ) = G(r 3 ,r 1 )u 2 P L (r 1 ,T 2 p'(r 2 ,r 3 ). (23) 



L Single- Degenerate version of Hb where the points ri and r 2 are merged, i.e. no 
scattering sequence takes place between these points: 

lL Singl S$) = 47rR2 J J ^i^ 2 ^ n (z2)C(^i)« 2 G 2 (ri,r 3 ) M 2 G ouf (r 2 ,R)Gl t (r 1 ,R) 

4:71 I I (z-[+Z2) ( z l+ z 2) 2 -l l \ 1 I \ 

= -p dz 1 dr 2 e- L ^e-^^G {r u r 3 )e lk " Ar2 - ri) e lk ™ t{ri - r2 \ (24) 

To analyze the respective contributions, we first investigate the dependence on 
the backscattering angle *& and, second, we perform angular integration to obtain 
information about the total backscattered flux j tot = J" dc ° s ^ 7(^). 

To evaluate the contributions from the respective processes we use a numerical 
Monte-Carlo algorithm simulating the propagation of the wave through the medium as 
a random walk. A comparable method has been used for example in [H] to analyze 
the backscattering cone for different shapes of the scattering medium. In all average 
Green's functions occurring in the above expressions, we use the ladder approximation 
l s = £ for the scattering mean free path. (Weak localisation corrections to l s will be 
considered in chapter [4] below.) 

The results of the angle-resolved analysis are shown in figure [3] For the Ha- 
contribution, we reproduce the characteristic backscattering cone (figure (3^a)). The 
contributions of the Lp- and Lsmgie-type scattering processes also yield positive peaks 
centered around i? = 0, but they turn out to exhibit a much wider angular distribution 
than the coherent backscattering cone and to be several orders of magnitude smaller 
in height. The only scattering processes yielding a negative contribution -- and 
hence being the only candidates for possible mutual cancellations among the different 
scattering scenarios - are ife-type scatterings which lead to a small cutback for every 
backscattering angle i? (figure |3^b)). 

In order to find out wether these contributions indeed cancel, we investigate the 
total backscattered flux 7 ioi of the respective contributions for different values of the 
disorder parameter k£. From the results obtained numerically, we derive the following 
asymptotic scaling laws in the limit of large hi: 

ltotHA ^(°- 42 )^ ( 25 ) 

Mltot HB ) ^-(O.lT) 1 ^, (26) 

Re( 7ioiiF ) ^(Q-19) 1 ^, (27) 

Re (^ Sl „ 9 J^(0.24)^p. (28) 



Both in figure p^ and in eqs. (25 28), we give only the real parts of the respective 
contributions since the imaginary parts are cancelled when adding contribution from 
the corresponding conjugate diagram (e.g. 7# c = 1*h b )- We find the leading 
order contribution to coherent backscattering to originate from the .ff^-type scattering 
processes, whereas all the scattering processes including one additional scattering event 
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Figure 3. Angle resolved backscattered flux 7(1?), see eq. ^, for the scattering 
processes shown in figure!!] System parameters: M — 30 and b = 10. (a) Characteristic 
sharp coherent backscattering cone ^h a (#) , see eq. (flS|) . (b) The contributions of all 
the processes involving one additional scattering event (real parts of 7# B ($), 7l f ($) 
and 7L Si „ glc (^), see eqs. (20 22 24 1) are found to be several orders of magnitude 



smaller in height and wider in angular distribution. Hg is the only process to contribute 
negatively and widespread over all angles d whereas Lp- and Lg ing i e -type scatterings 
yield positive backscattering cones centered around $ = 0. 



scale in higher order of l/(k£). In order to gain further insight as to the origin of 
the leading order contribution we now additionally distinguish the H ^-contribution 
with respect to the number N of scattering events inside the medium: We first look 
at scenarios, where the photon undergoes a crossed process including a true multiple 
scattering sequence (Ha,n>2), and analyze separately the possibility of double scattering, 
where the photon leaves the medium after but two scatterings {Ha,n=2) as shown in the 
top left sketch of figure [4] Thereby, we obtain the following asymptotic behaviour for 
large kt. 

ltot HAtN ^ -> (0-42)1, (29) 



Itoti 



->• (0.62) 



kt 
In {kt) 



(30) 



Hence, the double scattering case of the crossed sequence gives the leading order 
contribution to the backscattered flux [12]. This leading order term of the order 
l/{kt) can thus not be compensated by any scattering scenario including an additional 
scattering event, since they were all found to scale in higher order of l/(k£). 
Furthermore, we checked for cancellations among the contributions scaling like 
\n(k£)/(k£) 2 : 

,ln(H) 



itot„ A , N>2 + itat„ B + itot Hc -* Q-28 y 2 7 ^ 0, 



Jtot 



H A ,N>2 



+ 2Re(2 ltotHB + ltotLF + 



Jtot, 



Single 



\n(k£) , 
^0.80^A^ ^0, 



(31) 
(32) 



where (32) amounts to the sum of all processes (H B 



(k£) 2 
L D , L F and L Single ) depicted 
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in figure [l] and their complex conjugates. (Only Ha is identical to its own complex 
conjugate after exchanging r x and r 2 .) In both cases, it is not possible to establish 



any cancellation of the i?A,Af>2-type contribution. Since (31) corresponds to the the 



mechanism of cancellation proposed in [8], we state here a clear contradiction to the 



theory stated beforehand. Eq. (32) indicates, however, that including all other scattering 
scenarios with similar scaling behaviour is not the right track towards recovering flux 
conservation in coherent backscattering, either. 

4. Coherent backscattering and weak localisation 

As explained before, propagation processes leading to exact phase cancellation between 
the wave and its complex conjugate counterpart survive disorder average. In an infinite 
medium, loop propagation processes (i.e. scattering sequences with equal start and end 
point) leading to this very scenario are known to be responsible for the effect of weak 
localisation in the case of weak disorder: Just as for the ladder and crossed diagrams 
yielding the background distribution and the coherent backscattering cone, respectively, 
for a loop the phase shift between the wave and its complex conjugate vanishes both for 
the case of equal and reversed pathways. Including these interference paths between the 
wave and its complex conjugate in the description was found to influence propagation 
properties since it changes the scattering and the transport mean free path [IB] . 

The top right sketch of figure lil (H AN=2 ) shows the weak localisation scattering 
scenario including a crossed propagation process (Ha) for the case of double scattering 
(N = 2). However, if we treat weak localisation in a finite rather than in an infinite 
medium, now taking into account the presence of the boundary surface, we must allow 
the possibility for the following processes to occur: Starting and endpoint of the loop 
(located at points ri and r 2 ) may also be located on the other side of the boundary, 
thus outside the medium (figure 4 bottom left: H AN=2out and bottom right H AN=2in , 
respectively). In the first case, we thus obtain a further possibility for the wave to leave 
the medium and propagate to the detector (H AN=2out ), and in the second case a new 
way for the wave to enter the medium (H AN=2in ). Of course it is also possible for 
both, starting and endpoint, to find themselves outside of the scattering medium. By 
drawing the scattering processes corresponding to this very last scenario, we retrieve 
the diagram giving rise to the leading coherent backscattering contribution (figure [4] top 
left: H AN=2 ). 

Consistency now requires that, as soon as we account for one of the before mentioned 
propagation processes, we immediately have to consider the other possible constellations 
as well. This directly leads to a joint treatment of weak localisation and coherent 
backscattering in the frame of a full treatment of each possible scattering process in the 
presence of a boundary surface of a finite scattering medium. We will therefore now turn 
our attention to the analysis of this completed description of coherent backscattering 
which, as we will see, turns out to be crucial to ensure flux conservation in a finite 
medium. For this purpose, we will restrict ourselves to H AN=2 -type scattering in order 
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Ha, (n=2) 




A, (N=2), out 




Figure 4. Scattering processes taken into account in the frame of a complete 
description of coherent backscattering (here restricted to double scattering) together 
with weak localisation. These processes arise naturally when accounting explicitly for 
the boundary surface of a finite scattering medium: If the starting point r s and the 
endpoint rf of the propagation are both located inside the medium we deal with a 
weak localisation loop (top right: H AN=2 ), whereas the constellation where the end 
point is located outside the medium while the starting point is still inside (bottom left: 
) contributes to the propagation from r s to the detector. Similarly, if the 



I 

A,N=2.out 



endpoint is inside but the starting point outside (bottom right: H A N=2 in ), we obtain a 
new possibility for the wave to enter and propagate to 17 . The coherent backscattering 
process (top left: Ha,n=2) corresponds to the case were both, the starting and the 
endpoint, are located outside the scattering medium. The four processes cancel each 
other mutually, thus ensuring flux conservation. 



to establish a full description of the process yielding the leading contribution to coherent 
backscattering for weak disorder, as shown in chapter [3j 



5. Flux conservation for leading order processes 

The top right diagram of figure [4] gives rise to propagation from r s to rf with intermediate 
scatterers located at points ri and r 2 : 



P 



Ha, 



r s,r f; 



u 




dY\dY2G s iG S 2\Gi2\ Gr^G^/, 



(33) 



with G\2 — C(l r i — r 2|)? G s \ = G(\r s — ri|), etc., the averaged Green's function 
as given in Q. As shown in [13], this process is associated with a change of the 
scattering mean free path l s such that, in an infinite medium, the sum of Pjj A N _ 2 {r s , r f) 



and the single-step ladder propagator P ^(r s , rf), see eq. (15), is normalized, i.e. 
f drf(P 0L (r s ,rf) + Pjj a N - 2 ( r s, r f)) — 1> see a ls° ec i- (45) below. Similarly, the diagrams 



rl 



H Ain and H Aout (bottom right and bottom left, respectively), describe the following 
contributions to the incoming and outgoing wave intensity: 



Ph ( r f) 



P H ( r i 



U 



U 



rfr 1 dr 2 ^ n (z 1 )V'*„(^)|G' 12 | 2 G; / G 2/ (34) 

rfr 1 rfr 2 G sl G: 2 |G 12 | 2 G out (r 2 ,R)Gl t (r 1 ,R) (35) 
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with ip in and G out given by eqs. ( 12|l3 ). Connecting these building blocks and the single- 
step ladder propagator with each other gives rise to the following integral equation for 
the average wave intensity J(r) = |?/>(r)| 2 : 

/(r) = \^ n (z)\ 2 +P l H A , N ^Jr)+Jdr> (P 0ji (r,r') + ^ AiiV=2 (r,r')) I(r')(36) 

from which the outgoing intensity flux results as: 

7(0) = 7#A,iV=2 W + 7£*tended($) (37) 

where the first term represents the double scattering contribution to coherent 
backscattering, whereas the second term amounts to a renormalized background: 

lExtendedW = AnR 2 u 2 f dz I(r) ( \G mt (r , R) | 2 + P l H A> N=2 out (r) ) (38) 

In order to verify that the sum of both terms, eq. (37), fulfills flux conservation, we 
proceed as follows: For k£ ^> 1, it can be shown that the main contribution to the 
double integral over r 1 and r 2 in eq. (33) originates from the cases where ri and r 2 are 
located very close to each other. We may hence expand the phase-sensitive exponents 
in the Green's functions occurring in eq. (33) in first order of the difference e = r 2 — ri 
and obtain: 

e~ £ /£ 



P° HAN ^ 2 (r s ,r f ) ~ J drdeP , L (r s ,r)P , L (r,r { )u 2 ^ 



(4< 



2 



rfrPo,L(r s ,r)P , L (r,r f )^^ (39) 

where r = (r x + r 2 )/2 and q = k ( ?J _*, — Z a Zl\ ) ■ Its absolute value g = |q| = 

ky/2 — 2cosf? depends on the angle $ = <(r — r s , r — r f ). Eq. (39 ) can be interpreted as 
a two-step ladder propagation (first from r s to r, then from r to rf) with non-isotropic 
intermediate scattering event defined by the angular distribution: 

„. n . arctan (k£y/2 — 2 cos $) , , , 

P(cos^) = \ V ] . (40) 

kiy 2 — zcosv 

Similarly, for k£ S> 1, the contributions to the incoming and outgoing intensity, 

eqs. ( 34p5 ), can be written as combination of two ladder steps with angle-dependent 

scattering at r: 

4 A ,^>) * |^ m (.)| 2 P 0>L (r,r f )P (i^r) (41) 



Ph (n)^ / drP £ (r t , r) iGoutir, R) | 2 P ( ^ ^, ' ^ — *- 

The same approximation applied to the maximally crossed double scattering 
contribution yields: 

r L _ 

1Han .M ^ 4vrP 2 / dz |^(z)|V |G ou4 (r,R)| 2 P(cos^) 
Jo 

f ' OS ' ; -P(costf) (43) 



cos $ 
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\T2 




Figure 5. Anisotropic scattering sequence including ladder-like and i?A,Af=2-processes 
(see figure |4J) in the approximate two-step description, see eqs. ( 39p3 ). The 
distribution of the angles ■& is given by (44 1. The sum of all these sequences reproduces 
the leading-order contribution to the coherent backscattering cone jh a n=2 ($) together 
with the correspondingly renormalized background, see figure [6J 



for L ^> £. Integrating eq. (43) over all backscattering angles, we recover our previous 
numerical result, eq. (29 ), for the asymptotic behaviour lim^-^oo k£ f Q 
7r( v / 2-sinh" 1 l)/4~0.42. 



2 iHa,N 



.M 



Using eqs. (39 43), the transport process described by eqs. ( 36p7 ) reduces to 
an anisotropic random walk of a classical particle, see figure [5] For such anisotropic 
scattering processes, flux conservation is naturally ensured. At each scattering event, 
the particle may either be scattered isotropically (as it is the case if only ladder diagrams 



are considered) or anisotropically according to the angular distribution, eq. (40). (More 
precisely, eqs. ( 36p7 ) do not account for sequences with two subsequent anisotropic 
scattering events which, however, are negligible for k£ ^> 1.) Adding both cases, the 
total angular distribution results as: 



G(costf) 



+ F(cos$) 



(44) 



where the scattering mean free path £ s turns out such that the total distribution is 
normalized, i.e. J^ dc ° s,& G(cos$) = 1. For k£ > 1, this results in 

£ 1 ^^_j L ^ n fHki)\ 

£ 



2k£ V W 2 ) 

what agrees with the result found in 



(45) 



Finally, figure 6 shows the solution of eq. (37) for k£ = 30 and b 



10, 

obtained by Monte-Carlo simulation of an anisotropic random walk as explained 
above. As expected, the renormalized background jExtended(^) is slightly reduced 
with respect to the ladder background such that, as we have checked, the sum 
of ^Extendedi^) and the double scattering cone 1h an=2 (^) fulfills flux conservation: 



1. 
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Figure 6. Comparison of the incoherent background distribution obtained when 
only ladder propagation is considered [jl adder ($)] with the case where the ladder 
approximation has been extended [7 'Extended^), see eq. (38)] by the Ha, N=2~^oop 
processes shown in figureHJ System parameters: k£ — 30 and b — 10. For the Ha,n=2~ 
scattering process we also show the corresponding backscattering peak YlH A n=2 {&)] ■ 
Whereas adding this ^h a jv=2"P ea k on top of pure ladder propagation ^Ladder would 
violate flux conservation, the latter is restored when replacing ^Ladder by the slightly 
reduced 7 Bxte 



tended • 



6. Conclusion and outlook 

Concerning the mechanism restoring energy conservation in the effect of coherent 
backscattering, we have obtained the following results: The leading order of coherent 
backscattering, found to arise from i^jv^-scattering (double scattering event of the 
crossed contribution), cannot be counterbalanced by scattering processes including one 
additional scatterer, since the contributions of these scattering processes where found 
to scale in higher order of I /(hi). We could not establish any cancellations between the 
higher order processes, either, which contradicts the previous theory presented in |Sj. 

On the other hand, it turned out to be crucial to develop a full description of 
coherent backscattering in the frame of the loop propagation processes as considered in 
the description of weak localisation in infinite disordered media |13j . Considering all 
possible constellations of a weak localisation loop with respect to the boundary surface of 
the scattering medium, one finds additional scattering contributions which equally have 
to be taken into account. For the leading order contribution we present an approximate 
description of the scattering process as an anisotropic random walk. Within the frame 
of this approximation, it is possible to obtain a flux-conserving description of coherent 
backscattering. For the remaining higher order processes depicted in figure [TJ we raise 
the hypothesis that each contribution to coherent backscattering is in the same way 
cancelled intrinsically in the frame of a complete description accounting for all possible 
constellations in the presence of a boundary surface. A possible ansatz in order to prove 
this hypothesis might be a generalization of what is known as the Ward identity [I 
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ensuring flux conservation for scattering processes in infinite disordered media, to the 
case of a non-translationally invariant scattering medium in presence of a boundary 
surface. 

From an aesthetical point of view, it is certainly satisfying to dispose of a theory 
of coherent backscattering where flux conservation is intrinsically built in. This is the 
main motivation of the present paper. Practical consequences concerning, e.g., the 
comparison between theoretical and experimental coherent backscattering cones in order 
to determine properties of the scattering medium, remain to be investigated. In this 
case, effects like internal reflections due to a refractive index mismatch [15] or the finite 
correlation length of the random medium - which have been neglected in the present 
paper - must also be taken into account. It will be interesting to see to what extent 
the flux-conserving scattering diagrams identified in the present paper then lead to 
predictions that differ from other approaches (e.g. [8]), and whether this will enable a 
more accurate determination of scattering properties such as the transport mean free 
path. 
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